THE BEST CONSTANT IN A FRACTIONAL HARDY 

INEQUALITY 



KRZYSZTOF BOGDAN AND BARTLOMIEJ DYDA 



Abstract. We prove an optimal Hardy inequality for the fractional Laplacian 
on the half-space. 



1. Main result and discussion 

Let < a < 2 and d = 1,2, The purpose of this note is to prove the following 

Hardy-type inequality in the half-space D = {x = (x\, . . . , Xd) &W l : Xd > 0}. 

Theorem 1. For every u £ C C (D), 

(1) \ iU ^_~ y ^! f dxdy > K d<a j D n\ X ) xfdx, 

where 

rid, a - 



r(«±M) a2 a 
and (fTJ) fails to hold for some u £ C C (D) if Kd,a is replaced by a bigger constant. 

Here B is the Euler beta function, and C C (D) denotes the class of all the con- 
tinuous functions u : R d — ► R with compact support in D. On the right-hand 
side of ([1]) we note the infinite measure x d ~ a dx, where Xd equals the distance of 
x = (xi, . . . , Xd) € D to the complement of D. Analogous Hardy inequalities, in- 
volving the distance to the complement of rather general domains, and arbitrary 
positive exponents of integrability of functions u, were proved with rough constants 
in [16] (see also [33l[TU|T7]). Thus the focus in Theorem[T]is on optimality of Kd, a - 
We note that Kd,i = and k^q > if a ^ 1 (see the proof of Lemma [2]). 

Theorem [1] may be viewed as an application of ideas of Ancona [Ij and Fitzsim- 
mons [19] , Indeed, consider the Dirichlet form £ , with domain Dom(£), and the 
generator £, with domain Dom(C), of a symmetric Markov process ([22], [35], [29]). 
and a function w > 0, and a measure v > on the state space. The following result 
was proved by Fitzsimmons in 19J. 

(3) If Cw < — wu then £(u,u)> J u 2 dv , u S Dom(£) . 

Thus, every superharmonic function w (i.e. w > such that Cw < 0) yields a 
Hardy-type inequality with integral weight v = —Cw/w. For instance, in the proof 
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2 K. BOGDAN AND B. DYDA 

of Theorem[T]we will use w(x) — x d a ^ . Full details of ([3|), and a converse result 
are given in [T^J Theorem 1.9]. Recall that 

£(u, v) — ~{Lu, v) , if u G Dom(C) , v £ Dom{£) , 

([22], [35]). Therefore, equality holds in © if u = w G Dom(C), see OH (1.13.c)], 
j9}. If w ^ Dom(C), or £u> does not belong to the underlying L 2 space, then, as 
we shall see, the optimality of v = —Cw/w critically depends on the choice of w. 
According to |7J , the Dirichlet form of the censored a-stable process in D is 

c(«,»)=u_ tt / / m^m^i^m dxdy) 

2 JdJd \x-y\ a + a 

with core C^°(D) (smooth functions in Cd(D)), the Lebesgue measure as the refer- 
ence measure, and the following regional fractional Laplacian on D as the generator 

(0 (3.i2)], mimy- 

J Dn{ly _ x>£} \x-y\ d+a 
Here A d ,- a = T((d + a)/2)/(2~ a n d ^\T(-a/2)\), Clearly, {TJ is equivalent to 

(4) C(u,u)>Ad- a K d , a / u 2 (x) x~ a dx. 

Jd 

Recall that the Dirichlet form of the stable process killed when leaving D is 
K(u,v) = -A d ,- a / / dxdy, 

L Jm d jR d f ~ y\ 

with core C%°(D), the Lebesgue measure as the reference measure, and the frac- 
tional Laplacian (on Mr) as the generator, 

e->0+ 7R<in{|2/-a:|>£} F — 2/1 

(see, e.g., [7]). Decomposing R d = D U Z? c , one obtains 

fC(u, u) = C(u, u) + / u 2 (x)n D (x)dx , u€C™(D), 
Jd 

where (the density of the killing measure for D is) 

d-1 . 



k d (x)= A d ,- a \x-y\ a dy=-A d ,- a ^, a+ /, % 
Jd- n rfa±a\ 



« d '- a r(s±M) 



see (2.3), (5.4)-(5.6)]. It follows from gj) and Theorem Q] that 

i^r(^). r 



C(w,u) > .4d,-a(«d,a + TVa+rf , — )/ u (x)x d a cia: 

(5) = 1 2 - / u 2 (x)x- a dx, 



ID 

for all u £ C%°(D), and the constant T 2 (^^)/n is the best possible. 

We like to note that in some respects, the censored stable process is a better ana- 
logue of the killed Brownian motion than the killed stable process is (see [14j [32] , 
and (37] [31]). We suggest the former as a possible setup for studying Dirichlet 
boundary value problems for non-local integro-differential operators and the cor- 
responding stochastic processes ([30], [38]) on subdomains of K d ([3]), beyond the 
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"convolutional" case of the whole of M. d ([2lJ [4]). In this connection, we refer to 
[25l l26l [24] for Green-type formulas for the censored process. 

The reader interested in fractional Hardy inequalities may consult [MJ H7J (Ml 
1141 1161 fTT) . In particular, fl} improves a part of the (one-dimensional) result given 
in [33j Theorem 2]. The fractional Hardy inequality on the whole of M rf is known 
as Hardy- Rellich inequality, and the best constant in this inequality was calculated 
in [2HI ES] (see also [4] for Pitt's inequality). As seen in [16], the asymptotics of the 
measure dist(x, D c )~ a dx agrees well with the homogeneity of the kernel \y — x\~ d ~ a 
in ([1]). Noteworthy, if a < 1 and D is a bounded Lipschitz domain, then the best 
constant in (QJ is zero f[l6]). 

We like to make a few further remarks. Theorem Q] and the results obtained 
to date for Laplacian and fractional Laplacian suggest possible strengthenings to 
weights with additional terms of lower-order boundary asymptotics ([HI [23] [4]), 
and extensions to other specific or more general domains ([MIS])- To discuss the 
latter problem, we consider open f2 C D, and its Hardy constant, k(H), defined as 
the largest number such that 

Iff (u(x) — u(y)) 2 , , f u 2 (x) , „ 

o / / id+a dxdy> K (n) / -T^ dx, ueC c {n). 

2 JnJn \x~y\ d+a J n dist(x, n c ) a 

Note that k(O) > if f2 is a bounded Lipschitz domain and a > 1 ([16 ). Let 
u € C c (n) C C C (D). We have 

1 I I (u(x) - «(,))» < 1 f /■ («(«) - u(y)f 

and 



m 2 (x) , f u 2 (x) , 
w -dx> —^-dx , 



dist(a;,fi c ) Q ~ 7^ ^ 
thus k(51) < Kd.Q. We conjecture that k(£1) = a for a S (1, 2) and convex f2, see 
[51)1 Theorem 11] for case of the Dirichlet of Laplacian. 

Examining ([2]) we see that Hd.a —> oo if a — > 2. This corresponds to the fact that 
the only function u £ C c (Cl) for which the left hand side of ([1]) is finite for a = 2 is 
the zero function, see [T5][rB]. However, Ad- a ^d,a — ► 1/4 and r 2 (i^)/7r — > 1/4 as 
a — > 2, an agreement with the classical Hardy inequality for Laplacian ([H]) related 
to the fact that for u G C^{D), A^ 2 u — > Am and C{u,u) — > — J Au(x)u(x)dx = 
J | Vu(x)\ 2 dx as a — > 2 (the latter holds by Taylor's expansion of order 2, and a 
similar result is valid for K,). For the vast literature concerning optimal weights 
and constants in the classical Hardy inequalities, and their applications we refer to 

[a na nu eqi Ea mi e] . 

Our primary motivation to study Hardy inequalities for non-local Dirichlet forms 
stems from the fact that the converse of (J3j stated in [19l Theorem 1.9] allows for a 
construction of superharmonic functions, or barriers ([1]), when a Hardy inequality 
is given. These functions may then be used to investigate transience and boundary 
behavior of the underlying Markov processes (pQ, [7J, [TjJ, [H]). In particular, we 
expect that the results of [TH [T7J may be used to obtain, for the anisotropic stable 
([TD1E3) censored processes, the ruin probabilities generalizing [7J Theorem 5.10], 
and to develop the boundary potential theory on Lipschitz domains ( [Jj [26l [24] ) in 
analogy with those of the killed stable processes ( [13l [6] ) . We also like to mention 
the connection of optimal Hardy inequalities with critical Schrodinger perturbations 
and the so-called ground state representation [2T] , 
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Despite the general context mentioned above, the paper is essentially self- 
contained and purely analytic. In particular we directly derive Fitzsimmons' ratio 
measure by a simple manipulation with quadratic expressions, ©, not unrelated 
to the ground state representation of [21j (4.2)]. Theorem [1] is proved below in this 
section. In Section [2] we calculate auxiliary integrals. 

In what follows, \x\ = {x\ + • ■ ■ + x 2 ) 1 / 2 denotes the Euclidean norm of x = 
(xi, . . . , Xd) G R d , and B(x, r) denotes the Euclidean ball of radius r > centered 
at x. For d > 2 we occasionally write x = (x', Xd), where x' = (x\, . . . , Xd-i), and 
we let ||x'|| = maxt = i r . i( j_i \xk\, the supremum norm on R d_1 . 

Proof of Theorem^ For u,v G C^°(D) we define (Dirichlet form) 
and (its generator) 

r t \ v f u ^ ~ j 
Lu(x) = hm / -: rj- — ay, 

^+J Dn{ly _ x>£} \x-y\ d + a 

so that Ad,- a £ = C, Ad,- a £ = ^r) 2 \ and £(u,u) equals the left-hand side of ((T|). 
Let p £ (—1, a), x — (xi, . . . , Xd) G D, 

w p (x) = x p d . 

By (5.4) and (5.5)], 



a — i 

7T— r( 



(6) Cw p (x) = 7(a,p) (a —, — x d a w p (x), 

1 \ 2 > 

where the (absolutely convergent) integral 

( 7 ) 7{a,P)= n > 



o 



(1 -t) 1+c 

is negative if p(a — p — 1) > 0. Guided by the discussion in Section [1] we let 



(8) "0»0 = —fa- = -7(«,P) rW - d 

Since, for each t G (0, 1), the function 

(t p — 1)(1 — t a ~ p ^ 1 ) 



(1-t) 



l+a 



is convex and symmetric with respect to (a — l)/2, therefore p i— > 7(a, p) has a 
non-positive minimum at p = (a — l)/2. By Lemma [5] below, ([5]l. and Q, we 
obtain (JTJ for u G C^°(D) C Dom(C), with given by ([2]). The case of general 
u G C C (D) is obtained by an approximation. 

Since the setups of [19] and [7] are rather complex, we like to give the following 
elementary proof of (P). Let w = iU( a _i)/2, w G C C {D), x,y G D. We have 

/ / \ / «2 2/ ,w(y)-w(x) 2 tu(a;)-«;(y) 

(u(a;)-u(y)) + u (x) — + u (y) — 

w(x) w{y) 

(9) = w(x)w(y)[u(x)/w(x) — u(y)/w(y)] 2 > . 
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We integrate ((9]) against the symmetric measure l\ y ^ x \ >e \x — y\ d a dxdy, and we 
let e — > + . According to the calculations above, 

1 f f (u(x)~u(y)) 2 f 2 y w(x) - w(y) da; 

-D-D D {y€-D: a=|>e} 

D 

To complete the proof we will verify the optimality of Kd,a- In what follows we 
denote p = If a > 1 then we consider functions n„ such that 

(i) u n = lon [-n\n 2 ] d - 1 x [1, 1], 

(ii) suppu n C [-n 2 -l,n 2 + l] d - 1 x [-L 2], 

(iii) < u n < 1, |Vu n (x)| < exj 1 and |V 2 u„(a;)| < ex" 2 for x e D. 
If a < 1 then we stipulate 

(i') !)„ = 1 on [— n 2 , n 2 } d ^ 1 x [l,n], 
(ii') suppiv, C [— n 2 — n, n 2 + x [|,2n], 

(iii) < v n < 1, |V«„(x)| < cx^ 1 and |V 2 u„(a;)| < cx^ 2 for x e D. 
We define (for any a G (0, 2)), 

(10) itn(x) = u„(a;)a;^ ■ 

We have 

^ML dx > f ^dx^(2n 2 ) d - 1 logn. 



{x: ||x'||<n 3 , ±<x d <l} 

Thus, by Lemma[4] below, K ( i. a may not be replaced in ([1]) by a bigger constant. □ 



2. Appendix 



Lemma 2. For < a < 2, 



(11) 7(a,^) = -- 

2 a 



Proof. Since 
we are led to considering 



1 t P _ t a-l _ 1 + t a- p -l 



7(a,P)= / ^ _,u+ Q 



B„(a,6) = / ^(l-i) 6-1 ^. 
Jo 

Here and below a > 0, b > —2, and < k < 1. We will also assume that 6 ^ 0, 1. 
Using t a_1 = t a_1 (l — t) + t a , and integration by parts, we get 

B K (a, b) = ^B K (a, b + 1) - i K a (l - K ) b , 

hence 

a + b ( a+b + 1 „ , , , „ N 1 „,„ Nb+1 \ 1 6 



B«(o, &) = ^- [ b+1 B K (a, b + 2)- ^-^(1 - «0" +1 J - - k) 
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Clearly, 7(0, p) = lim^!- [B K (p + 1, -a) - B K (a, -a) - B K (1, -a) + B K (a - p, -a)] 
For a/lwe have, 

B K {p + 1, -a) - B K (a, ~a) - B K (1, -a) + B K (a- p, -a) = 1 x 

a{a — 1) 

{(p+ 1 - a)(p + 1 - a + l)B K (p + 1, 2 - a) - (a - a) (a - a + 1)-B re (a, 2 - a) 
-(1 - a)(l - a + 1)B K (1,2 - a) + (a - p - a)(a-p- a + l)B K (a-p, 2 - a)} 

+ a r_/p + j _ a wp+l + / Q _ Q wa + /j _ Q wl _ ( a _ _ a ) K a-P] 

a(o; — 1) 



All expressions in the square brackets, and their derivative, vanish at k = 1. Thus, 
they do not contribute to the limit as k — > 1. For a^lwc get 

1 



7(a,p) 



{(p + 1 - a)(p + 2- a)B{p + 1, 2 - a) 



a(a — 1) 

(12) -(1 - a)(2 - a)B(l,2 - a) + p(p - l)B(a - p, 2 - a)} . 

By the duplication formula T(2z) = {2-kY 1 ! 2 2 2z - 1 ' 2 T(z) T(z+l/2) with 2z = 2-a, 
for p = (a — l)/2, this equals 



3 — a „,a + 1 
-— 2-«) + l 



-r(^i)r(2-a)/r(^ ) + i 



-r(^-)r(^-)/r(-)2 1 -« + i 



5(^,^)2- -1 



2 ' 2 

We thus proved pip for a 7^ 1. The case of a = 1 is trivial. In fact, 7(1, 0) = 0. □ 
Lemma 3. Let — 1 < r < a < 2 and a > 0. There exists a constant c such that 

L js y i+« ^- m ~ a(avxd)r 

JD\B(x,a) \ x V\ 

for every a > and x £ D. 

Proof. Let a, t) = -B(x, f) \ -B(ir, s). If a > x d /2 then 



.'/,'/ 



£»\B(x,o) 1 35 ~ 2/1 



y' d 



k^O^ l £ | nB(a;,2 fc a : 2 fc + 1 a) (2 Cl) +a 
00 

k=0 



dy 



If a < Xd/2 then 



y' d 



I DC\B(x,a,x d ) \ X V\° 

and, by first part of the proof, 



— dy<cx r d a a , 



J d ld+a dy < cx r - a . 



D\B(x,x d ) 



This ends the proof. 



□ 



Recall that p = and u„ is defined by (TIT))) . 
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Lemma 4. There exists a constant c independent of n, such that 

: rj— dydx <cn ( > + 2n d , a \ u n (x) x d dx. 

djd \ x — y\ J 

D 

Proof. To simplify the notation we let K n = suppu„ and u = u n , v — v n . By @ 
and (JU) we have 

(u(x) — u(y)) 2 , , „ f 9 , . „ , 

\ x -y\d+a dX dy = 2Kd - a J ^ X d dX 
D D D 

(v(x) — v(y)) 2 , . . . , , 
' ■w(x)w(y)dxdy. 



D D 



\x-y\ M 



We will estimate the latter (double) integral by cn 2 ^ 1 ', by splitting it into the 
sum of the following six integrals I\ + . . . + Iq . 
We will first consider the case of a > 1 . 

If x <E K n and y G B(x, then \v(x) — v(y)\ < c\x — ylx^ 1 , as follows from 
(ii) and (iii). We thus have 



< 2 



> Jl'J w{x) W (y)dydx 

k„ Jb(x,^) f — y\ 

2p-2 



< c'n 2 < d - l >. 
A similar argument gives 

h = 



{V( : X) -^1 )2 W (x) W (y)dydx<cn^. 

{x:x„>i}iiJ(x,i) \x-y\ d+a 

We then have by Lemma [3] for a = 1/4 and r = p 

= / A . 

JD JD\B{x±) F 2/1 



< 



ld+^ w(x)w(y) dydx 



I K n J D\B{x,\) \x~y\° 

If d > 2 then we consider P n = {x e M d : ||x'|| > n 2 - 1 , < < |} and 
P° =F„n{ieK J : ||z'|| < n 2 + §}. We obtain 

(v(x) — v(y)) 2 , . , . , , 
— r£^- w(x)w(y) dy dx 
\x- y \ d+a y ) w> » 



h = 



— / / 1 C rrr dx 

ipo Jd\b(x,±) F - 2/1 

< c'|P> Q < cV^- 1 '. 
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We let R n = {x e R d : \\x'\\ < n 2 - 1 , < x d < \ } if d > 2, and we let 
Rn = {x e R : < x < 1} if d = 1. We have 



h = 



< 



(v(x) — v(y)) 2 , . , s , , 
I w{x)w{y)dydx 

lR n JD\B{x,^) F 2/1 



In the last inequality above we have used Lemma [3] with a = and r = p. 

We define L„ = {a; e R : ^ < s < |} in dimension d = 1, and for d > 2 we let 
L„ = {iel i : < n 2 - 1 , I < x d < ±}. We have 

(v(x)-v(y)) 2 



, , , 1N | T 7; id +Q ™{x)w{y)dydx 

I L n J DnB{x,^,\) \ x ~y\ 

w(x)w(y) 

lL„ J{y.0<y d <±} \ x - V\ C 

For d > 2 and x £ L n we have 

rfy /" 

(|x'-y'| 2 +2:2)(rf+«)/2 

{y:0< y(i <i} {y:0<y d <-i-} 



f dy <c f 

J \x-y\ d + a ~ C J 



71 \J {y'£M. d - 1 :\x'-y'\<x d } 
-a-1 



+ f "| *L 

J W eR^:\x'-y'\>x d} J iW - yV + X 2 )(1+*)/ 2 



< c 



n 

thus 



7 6 < c / ffir^H^^cv^- 1 ) 

Jl„ ^ n j n 



IL 

The case of d = 1 is left to the reader. 

We now consider the case of a < 1. We have 

(v(x) — v(y)) 2 , . , . , , 

^— p — -j^- w(x)w(y) dx dy 

\x-y\ d+a w w y 



7 = 



D 15 

< 



DJB(x,\) J{x:x d >%} JB(x,%) JdJd\B(x,%) J P n J DnB(i,if) 



'{x:0<:cd<2} JDnB(i,|,f) JL„ JflnB(x,i,f) 

= /l +/ 2 +/3 +Is +ie, 

where 

Pn = {xeR d : Hi'll >7i 2 -7i, 0<x (i < |}; 

L„ = {xeR d :\\x'\\<n 2 -n,2<x d <^}, 

for d > 2, and P n = 0, L„ = (2, 5) for d = 1. We estimate the integrals Jfc in a 
similar way as for a > 1. The details are left to the reader. □ 
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Similar but simpler estimates were given in [16j to prove that the Hardy constant 
of a bounded Lipschitz domain (e.g. of an interval) is zero if a < 1. We also 
like to mention that there is an alternative proof of Lemma U (not given here) , 
which explicitly uses the fact that that w 2 (x) = x^' 1 is harmonic ([7]) for A^ 2 . 
Similarly, the best constant, 1/4, in the classical Hardy inequality for the half-space 
D is obtained by considering w(x) = y/xd in Fitzsimmons' ratio v = —Aw/w. 

Acknowledgement. We thank Grzegorz Karch for a discussion on recent results 
on non-local operators, and Rodrigo Bahuelos, Malgorzata Kuchta and Jacek 
Zienkiewicz for discussions on classical Hardy inequality. 
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